Adiabatic preparation without Quantum Phase Transitions 
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Many physically interesting models show a quantum phase transition when a single parameter is 
varied through a critical point, where the ground state and the first excited state become degenerate. 
When this parameter appears as a coupling constant, these models can be understood as straight- 
line interpolations between different Hamiltonians Hi and Hy- For finite-size realizations however, 
there will usually be a finite energy gap between ground and first excited state. By slowly changing 
the coupling constant through the point with the minimum energy gap one thereby has an adiabatic 
algorithm that prepares the ground state of i/p from the ground state of Hi. The adiabatic theorem 
implies that in order to obtain a good preparation fidelity the runtime r should scale with the 
inverse energy gap and thereby also with the system size. In addition, for open quantum systems 
not only non-adiabatic but also thermal excitations are likely to occur. It is shown that - using only 
local Hamiltonians - for the Id quantum Ising model and the cluster model in a transverse field the 
conventional straight line path can be replaced by a series of straight-line interpolations, along which 
the fundamental energy gap is always greater than a constant independent on the system size. The 
results are of interest for adiabatic quantum computation since strong similarities between adiabatic 
quantum algorithms and quantum phase transitions exist. 

PACS numbers: 64.70.Tg, 03.67.Bg, 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

The promises of quantum computers in solving cer- 
tain difficult problems such as number factoring [l| and 
database search [2] have initiated a lot of research 
As no quantum system can be isolated perfectly from 
the environment the influence of decoherence is of 
great interest. In the conventional scheme of quantum 
computation, such errors have to be taken into account 
using expensive error correction schemes or by finding 
extremely fast quantum gates 0]. This has drawn at- 
tention towards alternative computation schemes such as 
measurement-based [1, 0] , holonomic 7] or adiabatic Q 
quantum computation that either circumvent error-prone 
parts of the computation process or even provide some 
intrinsic robustness against imperfect implementation or 
decoherence. 

The basic idea of adiabatic quantum computation 
(AQC) is to encode the solution to a difficult problem 
into the ground state of a problem Hamiltonian. A suit- 
able quantum system is prepared in the easily accessible 
ground state of a different initial Hamiltonian, which is 
then slowly deformed into the problem Hamiltonian. For 
slow evolutions the system will remain close to the instan- 
taneous ground state and will therefore finally encode the 
solution with high fidelity. The adiabatic theorem relates 
the maximum evolution speed with the spectral proper- 
ties of the time-dependent Hamiltonian Adiabatic 
quantum computation is polynomially equivalent to the 
conventional scheme of quantum computation and 
it is also believed to be somewhat robust to decoher- 
ent interactions at low reservoir temperatures pd| : In- 
tuitively, the ground state cannot decay and excitations 
can only occur for high temperatures. This is also backed 
by more formal calculations: For a time-independent sys- 



tem Hamiltonian and a thermalized reservoir with inverse 
temperature /3, the thermalized system density matrix 
with exactly the same inverse temperature can (in Born- 
Markov and secular approximation) be shown to be a 
stationary state of the master equation for the reduced 
system density matrix If for system Hamiltonians 
with a slow time dependence these approximations are 
also valid (Til . this would imply robustness of the 
adiabatic scheme as long as the reservoir temperature is 
smaller than the fundamental energy gap of the system. 

Adiabatic quantum algorithms bear strong similari- 
ties to the dynamics through a quantum phase transition 
[H,[I3)[IB|- They connect very different (macroscopically 
distinguishable) ground states and typically, the funda- 
mental energy gap is well controlled in the beginning and 
the end of the adiabatic algorithm. Somewhere in be- 
tween however, the fundamental gap may become very 
small and for nontrivial pro blems it will scale inversely 
with the system size Il7l|. Around the minimum en- 
ergy gap the ground state changes drastically ^8*1 and in 
the continuum limit n — > c», this would correspond to a 
quantum phase transition where an adiabatic evolution 
is impossible. However, as one will for AQC always be in- 
terested in a finite-size system, the degeneracy between 
ground state and first excited state is usually replaced 
by an avoided crossing. Here, the scaling of the energy 
gap with the system size is of great interest, as it will in 
practice mean a huge difference whether the gap scales 
inverse exponentially in n or merely polynomially. 

In this paper, two simple models will be studied 
that conventionally connect two Hamiltonians through a 
straight-line interpolation, along which in the continuum 
limit a quantum phase transition is encountered. It will 
be shown that with using a nonlinear path between the 
two Hamiltonians, the phase transition can be avoided. 
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II. ID TRANSVERSE ISING MODEL 

The one-dimensional quantum Ising model in a trans- 
verse field [l^ can be envisaged as a straight-line inter- 
polation 

H{s) = (1 - s)Hi + sHf (1) 

between the initial Hamiltonian Hi (at s = 0) and the 
final Hamiltonian (at s = 1) 



(2) 



where a^^^ denote the Pauli spin matrices acting on the 



■th 



qubit and periodic boundary conditions 



+1 



are assumed, which is however not essential for the scal- 
ing of the fundamental energy gap. Both Hamiltonians in 
([2]) obey a 180 degree rotational symmetry around all 
af-axes (bit- flip), which transforms (t| to —erf. 

The ground state of the final Hamiltonian is two-fold 

degenerate, with 

where cr^ |0) = -I- 



*[,^^) = I0...0) and 



(2)\ _ 



= 1. 



1) 



0) and CT^ |1) = — As the eigenval- 
ues of Hp are proportional to the number of kinks in 
the basis states trf \zi . . . z„) = (—1)^' \zi . . . Zn), the two 
ground states are well separated from the excited states. 

The unique ground state of the initial Hamiltonian in 
^ is given by 



1 



E 

2 = 



(3) 



where cr^ |— >) = +|^). The excited states (with at 
least one spin flipped) are again well separated from the 
ground state. 



Due to 



0, the bit-flip-parity is a con- 



Hf, erf 
i=l J 

served quantity during the interpolation ([T]), and since 
([3]) is even under bit flip, the flnal system state will for 
adiabatic evolution [l^, [20l . [2l| choose at s = 1 the even 
subspace (symmetry breaking) with the huge Schrodinger 
cat ground state = l/\/2 [|0 . . . 0) + |1 . . . 1)], if 

the system is initialized in ([3]) at s = 0. 

Finally, note that the overlap between initial and flnal 
ground state is exponentially small in the system size n. 



with the fermionic operators 7fc and single quasi-particle 
energies 



ek = ^/^-■icos^{ka/2)s{^-s), 



(5) 



see figure [T] left panel. The wavenumber ka covers the 
range ka = ^ {I + 2Z) : \ka\ < n. 

From equation ([5]) one obtains a minimum gap of 
.9min = 2sin (^) = 0{l/n} and indeed one also observes 
for constant speed interpolation s — t/r a quadratic scal- 
ing of the adiabatic runtime Tad — 0{n^} [l^ with the 
number of qubits n. At the critical point Scrit = 1/2, the 
Ising model undergoes a second order quantum phase 
transition jlSj . which manifests itself in a divergence of 
the second derivative of the ground state energy. 



B. Nonlinear Interpolation Path 

The linear interpolation scheme in equation ([1]) is 
rather simple but not necessary for the adiabatic theorem 
to hold [23, ll^] . The nonlinear interpolation we will con- 
sider here will connect the Hamiltonians of ([2l) by a series 
of straight-line interpolations along the Hamiltonians 



^^0 

Hi 



-J^at^Hi, 

n 



i=3 



Hk = -^(T^a-_^_i 



i=k+2 



- E ^i'^Ui ' 
i=l 
n 



(6) 



such that formally the time-dependent Hamiltonian is 
given by 



A. Straight-Line Interpolation 

At any point during the interpolation s € [0,1], the 
Hamiltonian H{s) can be diagonalized: The successive 
application of Jordan- Wigner, Fourier and Bogoliubov 
transformations [l^ H^l map the spin-l/2-Hamiltonian 
to a fermionic one, which (in the subspace of even bit-flip 
parity corresponding to even numbers of quasi-particles 
only) reads 



1 



(4) 



H{t) - EQ[^'=(*)]Q[i-^fc(*)] 



k=0 



{[l-Skit)]Hk + Sk{t)Hk+i} 



(7) 



where Q{x) denotes the Heavyside step func- 
tion and Sfe(t) = encodes a constant 
speed interpolation with nAt — t such that 
Sfe[(fc + l)At] = 1 = 1- Sfe+i[(fc + l)At] and 
H{t) = Hp. 

Given a spin chain with (constant) cr|^ cTj^^j^-interactions, 
one has to apply a correspondingly stronger external fleld 
to achieve an effective decoupling of the spins. Physically 
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the above nonlinear scheme could be approximated by a 
strong transverse magnetic field that within the distance 
between two spins rises linearly from zero to maximum 
and then travels at constant speed along the spin chain. 
Conversely, one might imagine to slowly pull the spin 
chain out of a region with a strong stationary transverse 
field. Recent proposals for quantum simulators of the 
Ising model in a transverse field 26] however, even 
provide a more direct control of the local spin-spin inter- 
actions, such that all regions of the phase diagram can in 
principle be reached. 

Note that the interpolation path ^ does not destroy 



diagonalized easily. For the Ising model we consider the 
controlled NOT-gate 

s^""""^ = ^ (i» + + l (1. - ® , (11) 

which is a Hermitian controlled- unitary operation (3], 
where i is the control and j the target qubit. It is 
straightforward to show that 



the bitflip symmetry, since 



Hk, erf 



e=i 



0. It is easy 



[a^ 1,] 



1.7, 



£cl cCNOT 



(12) 



to see that the ground state of a Hamiltonian in ^ 
(within the subspace of even bit-flip parity) is for 1 < 
fc < n — 1 given by 



1 



x/2 



[I+Of... fJ^+i] nk+2 . . . H„ |0 . . . 0) ,(8) '^fe+l,fe+2^fe(sjifc+l,fe+2 



where Tik = (cr^ + cr^ ) denotes the Hadamard gate 
on qubit fc, such that the overlap between two succes- 
sive ground states yields ^^'^ = l/\/2 for 
< fc < n — 2. In the last interpolation step, the ground 
state is even invariant ( = 1- Hence, 



-1 l^n / 

in every single step only slight transformations of the 
ground state are performed and intuitively, one may ex- 
pect that adiabatic preparation along this modified path 
should be more efficient than in the conventional scheme. 

For the first interpolation step Hq{s) = {1 — s)Hq + sHi 
in ([6]), the first two qubits evolve independently from the 
rest of the system and it is straightforward to obtain the 
eigenvalues from the nontrivial contribution of their four- 
dimensional subspace 



AS = -V5s2^^87+4- (n-2) + 2K, 

\1 = -s-(n-2) + 2K, 

A2 = +s - (n - 2) + 2k , 

= V5s2 _ 8s + 4 - (n - 2) + 2k , (9) 

where k G {0, 1, . . . , (n — 2)} counts the excitations re- 
sulting from qubits 3 . . . n, see figure [T] right panel. By 
continuity we identify Aq and A3 as belonging to the even 
subspace and thus one obtains in the first interpolation 
step for the minimum gap ff^™ = min [A3 — Ag] =4/ V5 
at Sinin =4/5. 

For the intermediate steps we have 



Hk{s) = {l-s)Hk + sHk 



k 



=1 



(1 



■i+l 



'k+2 



+ 1 



! = fe+3 



(10) 



in (l6|) for 1 < fc < 71 — 2. In order to obtain the spectrum 
conveniently, we will map the Hamiltonian using a suit- 
able unitary transformation to another one which can be 



Using the CNOT transformation at the transition region, 
the Hamiltonian (fTO|l is mapped to 



i=l i=k+3 
-(l--'^)'^fe+2-S^I+2(13) 



where it is visible that the qubit sets {l,...,(/c-fl)}, 
{{k + 2)}, and {{k + 3), . . . ,n} are mutually decoupled. 
Evidently, one obtains for qubits 1 ... (fc -I- 1) just the 
eigenvalues of the Ising model with open boundary condi- 
tions (i.e., with minimum energy — fc for a two- fold degen- 
erate ground state and fundamental energy gap 2) and for 
qubits (fc -|- 3), . . . , n a minimum energy of — (n — fc — 2) 
for the unique ground state and a fundamental energy 
gap of 2. The nontrivial part of the spectrum arises from 
the subspace of qubit (fc -I- 2), where one obtains for the 
eigenvalues 



A;^ = ±^1 - 2s(l - s) - (n - 2) 



2k. 



(14) 



where k g {0,1, . . . ,{n — 2)}, see figure [T] right 
panel. Therefore, the fundamental energy gap equates 
to gj^™ = V^, which is independent on the system size. 

The last step is trivial: As [7?„_i,if„] = 0, both 
Hamiltonians can be diagonalized simultaneously and the 
spectra of Hn-i and iJ„ are therefore just connected by 
straight lines as is also visible in figure [T] right panel. 

Therefore, we obtain with the new interpolation path 
a lower bound on the minimum fundamental energy gap 



even 



> %/2 



(15) 



in the even subspace. Note that as a sanity check, in 
figure [U right panel a perfect agreement is found between 
analytically and numerically calculated eigenvalues (27j 
is found. 

In the absence of a reservoir one would have an adia- 
batic scheme to efficiently prepare large Schrodinger cat 
states: If the quantum system is initialized in its initial 
ground state ([3]) (by applying a strong magnetic field) 
and then slowly evolved through the series of interpola- 
tions towards i/p, the adiabatic theorem will guarantee 
that the system will remain close to its instantaneous 




FIG. 1: [Color Online] Left: Lowest eigenvalues for the Id quantum Ising model in a transverse field ([2]) and linear interpolation 
U]) for n = 10. Thick lines show the analytical result for the eigenvalues of the even subspace ([SJ, whereas the thin dotted 
lines give the numerical result [S^l for all lowest eigenvalues. Right: Comparison of numerical calculation (dense symbols 
visible as thick lines) with the analytical prediction for n — (solid black lines, only shown in intervals [0, 1], [4, 5], [9, 10]) for 
stepwise interpolation ^ with n = 10 and r = 10. The different interpolation regimes are separated by dashed vertical lines. 
Throughout the interpolation scheme the minimum gap in the even subspace is always larger than y/2. The ARPACK package 
[27I ] only determines one part (e.g., the lower) of the spectrum, which is why not all numerical eigenvalues are shown. 



ground state. Evidently, the last step in ^ can also 
be omitted without changing the ground state, just the 
energies will differ. This docs not mean that one can 
adiabatically deform Hi in Hp in ([2]) in constant time, 
as the number of steps in ^ scales linearly with the 
system size n. In order to obtain a fixed final fidelity 
of the desired state |^'fi„ai) = l/\/2 [|0 . . . 0) + |1 . . . 1)], 
the adiabatic runtime Tad will scale nearly linearly with 
the system size Tad oc nlogn for the new scheme © 
instead of quadratically [l^ for the straight-line interpo- 
lation scheme ([T]). The logarithmic corrections are ex- 
pected from the increasing number of degeneracies in the 
first excited states Q]. However, this favorable scaling 
could also have been achieved by using the conventional 
scheme ( [1]) w ith straight-line but adaptive-speed inter- 
polation |28l . [29| . Therefore, the advantage of the non- 
linear scheme is rather given by the fact that constant- 
speed interpolation suffices and - more importantly - for 
couplings to a reservoir that do not destroy the bitfiip 
symmetry the new scheme provides some resistance to 
thermal excitations as long as fesTres < C*!!}. 



additional (time-independent) Hamiltonian of the form 



2 



e=i 



(16) 



to safely separate the even and odd subspaces. Unfortu- 
nately, it involves n-qubit interactions and is thus pre- 
sumably difficult to implement. In the following, we will 
therefore consider a model with a local Hamiltonian that 
has a unique ground state. 



III. PREPARATION OF THE CLUSTER STATE 

The cluster state on a lattice £ can be encoded in the 
ground state of the Hamiltonian 



(17) 



For realistic systems however it is to be expected that 
the conservation of bit-fiip parity will be destroyed by 
the existence of a reservoir, such that amplitude from the 
even s ubsp ace may leak into the subspace of odd bitfiip 
parity [UllOj (such effects may also limit the robustness 
of other models against decoherence, see e.g. [3l|). The 
net effect would be the decay of the desired Schrodinger 
cat state. In principle, this could be avoided with an 



where v ^ ^ denotes all sites v that are connected to 
the site ji in the lattice £ by a link. In order to see that 
the cluster state is the ground state of (fTT]) . we use the 
controlled Z-gate 

- \ (1. + ® 1, + ^ (1^ - ® = Sf , (18) 
which acts as 



(19) 



5 



in the computational basis. Obviously, it is Hermitian 
and [5**^^,5'^;^] — 0, such that we can unambiguously 
define the global unitary operation 



5. = n 



=4, 



(20) 



With using the relations (see also [32l |l 

S'^ [at ® Sgz = 1. « 
S'^ [at ® 1,] S?^ = at € 



(21) 



it is easy to see that the unitary transformation (|20p maps 
the Hamiltonian (fT7|) into decoupled spins 



ScHS 



(22) 



This implies that the ground state of (|17p is unique and 
separated by a safe energy gap of AE = 2 from the first 
excited states. From the ground state of (|^^ it can be 
obtained via 



2"-l 



^ ^(-1)- , (23) 



z=0 



where L^i, — 1 \i fi ^ v and L^^ = otherwise, and 
n = \C\ denotes the number of lattice sites. Above equa- 
tion already gives the recipe of preparing the cluster state 
using only the two-qubit unitary operations (|18p . How- 
ever, considering the destructive influence of an environ- 
ment it would be desirable to encode the cluster state in 
the ground state of a Hamiltonian, which could then be 
reached by adiabatic transformation, see below. 



A. The Id cluster state 

In analogy to the Ising model in a transverse field ([2]), 
one could prepare the cluster state on a Id spin chain 
adiabatically by a straight line interpolation ([T]) between 



Ht = 



and 



Hy = -t 



i=2 



(24) 



what is referred to as 'transverse field cluster model' in 
[3^ . It has been shown that the above model exhibits 
a quantum phase transition at Scrit = 1/2 [ssj, which 
can be understood as p4)) can be mapped via a duality 
transformation to two Ising models in a transverse field 
[3^ . Accordingly, this kind of preparation would be vul- 
nerable to nonadiabatic and thermal excitations in the 



large n limit, as the minimum energy gap will scale as 
0{l/n}, see also figure |3] left panel. 

Therefore, we will consider the stepwise interpolation 
scheme ([7]) along 



Hi 



1=1 



-a-^a^ 



n 

„z ^x „z „x \ „x 



i=4 



Z \ ^ Z Z ^Z 'X \ ^ 



i=k+2 



- E 



(25) 



The physical implementation of such a model would 
of course be more demanding: Given that the neces- 
sary three-body interactions can be simulated by logi- 
cal qubits (possibly built of several physical ones along 
the lines of |34|). the region of a non- vanishing transverse 
field gradient should here also be confined within the dis- 
tance between two logical qubits. In order to simplify the 
notation, we have written the intermediate Hamiltonian 
of (pS)) as in ^Tl} and just modify the number of existent 
links in the lattices Ck, see figure [51 Then, we obtain for 
















FIG. 2: [Color Online] Stepwise buildup of the cluster state 
on a Id spin chain with n = 7 spins within 6 steps using a 
runtime of r = 0{n}. In lattice £,1, all spins are unconnected 
(dark-blue) and in lattice £,2 only two spins are connected 
(grey-red) by solid lines. Thin dotted lines show the intended 
order of inclusion in the connected sublattice. Dark shaded 
regions (brown) show the presence of a strong transverse field, 
whereas in the light-shaded regions (yellow) the field is grad- 
ually decreased from maximum to zero. 

an intermediate interpolation step 

Hk{s) = (1 - s)Hk + sHk+i 



'\"2 



k 

^x z 



E< 

-(l-s)f7fefT^+l-(l-s)fT^+2 



''k^k+\^k^2 



^^k+\^k+2 



(26) 



6 



where it becomes visible that by applying the global uni- 
tary only on the existing links in lattice Ck, the model is 
again decoupled nearly completely 



E 



i=l : j^{(fe+l),(fe+2)} 



a 



k+2l 



(27) 



In above equation, only qubits {(fc+1), (fc+2)} contribute 
the nontrivial parts of the spectrum and we obtain for the 
eigenvalues 



A5 



~2^\ - 2s(l - s) - (n - 2) + 2k , 
\^ = -(n-2) + 2K, 

+2^/1 - 2s(l - s) - (n - 2) + 2k (28) 



with K e {0, 1, . . . , (n — 2)}. Note that the first interpo- 
lation step yields the same eigenvalues, as is also visible 
in the numerical calculation of the spectrum, see figure [3] 
right panel. 



The 2d cluster state 



In 2d, the cluster state can serve as a universal 
resource for measurement-based quantum computation 
_&]. Therefore, we will here also state some known 
facts about its adiabatic preparation via a straight-line 
interpolation ([T]) between the Hamiltonians 



f^ec 



(29) 



where C denotes a square lattice, such that five-body in- 
teractions are required in order to implement H^. For the 
direct straight-line interpolation one can map the above 
model with a duality transformation (3^ to one with a 
known quantum phase transition at Scrit = 1/2 [ssj . 

It is easy to see that one can choose an interpolation 
scheme where either only one or two links are added at 
once, see figure^l Whenever only a single link is inserted, 
we have the same situation as in the Id case, therefore 
we will consider the more interesting case of adding two 
links at once here. Suppose a node a is unconnected in 
lattice Ck (i.e., at s = with Hamiltonian Hk) and shall 
be connected to two other nodes (3 and 7 in lattice Ck+i 
(i.e., at s = 1 with Hamiltonian H^^i). Suppose also 
that the nodes (3 and 7 are already connected to other 
nodes in lattice Ck as depicted in figure ID Then, the 



corresponding Hamiltonian can be written as 

Hk{s) = {1 - s)Hk + sHk+i 



-{l-s) 



\ 



.(30) 



1/ ~ 7 



Using the unitary transformation Sc^ 1 we again arrive 
at a nearly decoupled model (where only three spins are 
mutually coupled) 



-s{alala; + alal + ala^^) . (31) 



Evaluating the eigenvalues of the 8x8 matrix generated 
by qubits {a, /?, 7} yields (for simplicity we only give the 
two lowest eigenvalues) 



5 - 10s(l - s) 

-{n - 3) + 2k, 
= -1 - (n- 3) + 2k, 



2s)2-K25s2(l-s)^ 
(32) 



where k e {0, 1, . . . , (n — 3)}. Therefore, we obtain a 
fundamental energy gap of ^min = \/5 — 1 at Smin = 1/2, 
which is independent on the system size. 

Note that we have so far allowed the external field to 
vary at a single site on the lattice only. If one does not 
have this constraint many of the steps can be performed 
in parallel, such that the overall time required to generate 
the cluster state is further reduced, compare also [l^l . 



IV. ADIABATIC QUANTUM COMPUTATION 

In the previous sections it was shown for specific mod- 
els that quantum phase transitions encountered along a 
straight-line path can be avoided if a nonlinear path is 
chosen. It would be great if this could be extended to 
adiabatic quantum computation, where for the nontriv- 
ial models numerical calculations based on a straight-line 
interpolation point to an inverse scaling of the fundamen- 
tal energy gap with the system size [ij, [13, [H, [13, [ll] . 
Not only would such a path usually admit a signifi- 
cantly shortened runtime of the adiabatic algorithm but 
also concerning robustness against thermal excitations 
a constant lower bound on the fundamental energy gap 
would be highly desirable. As a first step, we will show 
that there exist paths consisting of a series of straight- 
line interpolations, along which the fundamental gap 
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FIG. 3: [Color Online] Left: Numerical Calculation of the lowest eigenvalues for the adiabatic cluster model in a transverse field 
(|17|) for n = 10. Thick lines show the lowest five eigenvalues. Right :Numerical Calculation of the lowest eigenvalues for the 
stepwise cluster model (|25p for n — 10 and At = 1. Analytical predictions (solid lines) resulting from psp for k = perfectly 
match the numerical results. 






















FIG. 4: [Color Online] Stepwise buildup of the cluster state on 
a 2d spin lattice with n = 25 spins in 24 steps. Thin dotted 
lines show the intended order of inclusion in the connected 
sublattice. In each single step, either only one or two links 
are added. With regard to realization, one could imagine 
a region with a strong transverse magnetic field (dark-brown 
shaded region), which is turned off slowly (light-yellow shaded 
region) along a rectangular zig-zag trajectory as outlined by 
the thin dotted lines. 



can be calculated analytically. As one paradigmatic ex- 
ample for the problem class NP, we will illustrate this 
for Exact Cover 3 (ECS), which is a special case of 3- 
SAT. The EC3-problem can be defined as follows [|: 
Given m clauses each involving positions of three bits 
Ci = {pl,pf,p^) with 1 < i < m and 1 < pf < n one is 
looking for the n-bit bitstring {bi . . . 6„} with bi G {0, 1} 



that fulfills for each clause bpi 



■ bp2 + b„3 = 1 , where 



denotes the integer sum. A problem Hamiltonian encod- 
ing the solution to the EC3 problem in its ground state 
can be given as [H, [IJ 



Hp = 



1- L 



(33) 



where Zfe = i (1 - cr|). 

Unfortunately, since the neighborship topology in the 
irregular lattice defined by riij in an EC3 problem can- 
not be represented by next neighbor interactions for a 
fixed dimensionality and does not even have a regular 



structure, there is no obvious unique way of defining a 
nonlinear interpolation scheme, compare also [23l ]. Here, 
it will be proposed that by turning on the clauses one 
by one in p3|) whilst maintaining a unique ground state 
one should obtain a considerable improvement in the fun- 
damental energy gap for the average EC3 problem. We 
define a partial problem Hamiltonian via 



-Pi 



(34) 



The series of Hamiltonians in the nonlinear path can then 
be written as 



Hk = l- 



(35) 



connected by stepwise linear interpolations ([7]) for 
< A: < m and where the unique ground states are the 
symmetrized superpositions of solutions to a subset of 
clauses 



1*0 



1 



z ■ (ffi+...+H*)|z)=0 



(36) 



where we have denoted the ground state degeneracy of 
i/p + . . . + i/p by Nk . This of course does not uniquely 
define an interpolation scheme, since we may exchange 
the order of clauses without changing the final solution. 
For the linear interpolation between two projector Hamil- 
tonians Hk{s) = 1 - (1 - s) \^k) (^'fcl - s 
the simple result for the adiabatic Grover algorithm [38[ 
may be easily generalized [s^ to obtain the eigenvalues 



Ao 



/I = 2 ^^l±yi-4s(l-.sj 1^1- 

and \2 — . ■ . — A2"-i = 1, such that the min- 
imum fundamental energy gap equates to 

\{-^k\^k+i) 



1 



5i 



^Jrl (where < 1 
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is the reduction factor in the number of solutions at step 
k). For an average ECS problem with a unique solution, 
this would be a significant improvement, since typically 
a modest number of states will violate only a single 
clause. For extreme problems however, an exponential 
number of basis states could possibly just violate only a 
single clause, such that even re-ordering of the clauses 
would not yield a polynomial scaling of the fundamental 
energy gap. However, for direct interpolation between 
Hq and Hjn (where we would reproduce the adiabatic 
Grover algorithm [38| ) we obtain (for a problem with 
a unique solution such that N„i = 1) the minimum 
fundamental gap of g^l^'^^ — 1/V2", which is always 
significantly smaller than the smallest gap that could be 
encountered with stepwise interpolation. 

A further drawback of the scheme is that the Hamilto- 
nians l|35p are generally nonlocal and cannot be efficiently 
implemented, as can for example be seen from the formal 
representation of its ground state 



\S) . (37) 



Apart from the (generally unknown) normalization fac- 
tor A^fc, one would from inserting p4p in above equation 
obtain an exponential number of terms involving up to n- 
body interactions and also using three-body interactions 
in ([54|) as in the original approaches [1, [l^, [3] would not 
cure the problem. 

However, it is also known that the conventional 
straight-line adiabatic algorithms with local Hamiltoni- 
ans are more efficient than straight-line adiabatic algo- 
rithms that use a projection operator po| - even when 
both use the same initial and final ground states. It 
would therefore be an interesting option to approximate 
the projection operators in ([55)) by few-body Hamiltoni- 
ans with the same (or similar) ground states or to explore 
other nonlinear interpolation paths for local Hamiltoni- 
ans. 



we have shown that a nonlinear path provides a con- 
stant lower bound for the fundamental energy gap. For 
closed quantum systems, this would enable a very effi- 
cient adiabatic preparation of highly entangled ground 
states [l^ [2^. The used Hamiltonians in sections |TT] 
and IHII are always local and require local control of 
the couplings. They do not require complicated time- 
dependencies of the coupling constants as in adiabatic 
rotation [23| or with local adiabatic evolution [H, [ssj . 
For the quantum Ising model in a transverse field, the 
robustness of the scheme in the presence of a reservoir is 
hampered by the two-fold degeneracy of the final ground 
state - unless one is able to separate the even and odd 
subspaces by a safe energy barrier. For the cluster model 
in a transverse field one does not have this problem, since 
the ground state is always unique. Unfortunately, it re- 
quires 3-body interactions for the Id lattice and 5-body 
interactions for the 2d lattice. However, it is also known 
that the 2d cluster state can be approximated by the 
ground state of Hamiltonians with two-body interactions 
only [34J . The approach could also be of interest for other 
models with just two- body interactions having highly en- 
tangled ground states [U (4^ . Adiabatic quantum com- 
putation might also strongly benefit from using nonlinear 
paths in both algorithmic performance and robustness 
against thermal excitations. Beyond this, the presented 
idea may give rise to a class of Hamiltonians that can be 
used to calibrate numerical methods such as DMRG [431 . 
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